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Abstract 

We survey the recent development of the DeGiorgi-Nash-Moser-Aronson type theory for 
a class of symmetric jump processes (or equivalently, a class of symmetric integro-differential 
operators). We focus on the sharp two-sided estimates for the transition density functions (or 
heat kernels) of the processes, a priori Holder estimate and parabolic Harnack inequalities for 
their parabolic functions. In contrast to the second order elliptic differential operator case, the 
methods to establish these properties for symmetric integro-differential operators are mainly 
probabilistic. 

1 Introduction 

Second order elliptic differential operators and diffusion processes take up, respectively, an central 
place in the theory of partial differential equations (PDE) and the theory of probability. There are 
close relationships between these two subjects. For a large class of second order elliptic differen- 
tial operators L on M", there is a diffusion process X on M" associated with it so that £ is the 
infinitesimal generator of X, and vice versa. The connection between £ and X can also be seen as 
follows. The fundamental solution (also called heat kernel) for £ is the transition density function 
of X. For example, when = \ Y^j=\ ^ af") ' where (ajj(x))i<jj<„ is a measurable n x n 

matrix-valued function on M" that is uniformly elliptic and bounded, there is a symmetric diffu- 
sion X having £ as its L^-infinitesimal generator. The celebrated DeGiorgi-Nash-Moser-Aronson 
theory tells us that every bounded parabolic function of £ (or equivalently, of X) is locally Holder 
continuous and the parabolic Harnack inequality holds for non-negative parabolic functions of £. 
Moreover, £ has a jointly continuous heat kernel p(t, x, y) with respect to the Lebesgue measure on 
M" that enjoys the following Aronson's estimate: there are constants > 0, /c = 1, • • • , 4, so that 

ci p'^(t, C2|x — y|) < ]3(t, X, y) < C3p'^(t, C4|x — y|) for t > and x, y G M"". (1-1) 

Here 

V\t, r) := exp(-rVt). (1.2) 

See |36j for some history and a survey on this subject, where a mixture of analytic and probabilistic 
methods is presented. 
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Recently there are intense interests in studying discontinuous Markov processes, due to their 
importance both in theory and in apphcations. Many physical and economic systems should be 
and in fact have been successfully modeled by non-Gaussian jump processes; see for example, 
[71 [271 [28l [33] and the references therein. The infinitesimal generator of a discontinuous Markov 
process in M" is no longer a differential operator but rather a non-local (or, integro-differential) 
operator. For instance, the infinitesimal generator of an isotropically symmetric a-stable process in 

with a S (0, 2) is a fractional Laplacian operator c/^^^"^ := — c (— A)"/^. During the past several 
years there are also many interests from the theory of PDE (such as singular obstacle problems) to 
study non-local operators; see, for example, [1U[ and the references therein. 

In this paper, we survey recent development of the DeGiorgi-Nash-Moser-Aronson type theory 
for the following type of non-local (integro-differential) operators C on M": 

X JL^ Q / du(x) \ f 
'Cw(x) = - V — aij(x)— ] +lim {u{y) - u{x))J{x,y)dy, (1.3) 

^i^^9xi\ dXj J elO J{y^Kn,\y_^\y,y 

where either (ajj(x))i<jj<„ is identically zero or (aij(x))i<jj<„ is a measurable nxn matrix-valued 
measurable function on M" that is uniformly elliptic and bounded, and J is a measurable non- 
negative symmetric kernel satisfying certain conditions. Associated with such a non-local operator 
C is an M"-valued symmetric jump process X with jumping kernel J{x, y) and with possible diffusive 
components when (ajj (x))i<jj<„ is non-degenerate. Note that the jumping kernel J determines a 
Levy system for X, which describes the jumps of the process X: for any non-negative measurable 
function / on x M" x M", t > 0, x £ M" and stopping time T (with respect to the minimal 
admissible filtration of X), 



^fis,Xs-,Xs 



s<T 



^x 



f{s,Xs,y)J{Xs,y)dy ds 



(1.4) 



Our focus will be on sharp two-sided heat kernel estimates for £ (or, equivalently, transition density 
function estimates for X), as well as parabolic Harnack inequality and a priori joint Holder conti- 
nuity estimate for parabolic functions of C. When (ajj (x))i<jj<„ = and J{x,y) = c\x — 
for some a £ (0,2) in (jl.3p . £ is a fractional Laplacian ciA"/^ on M" and its associated process X 
is a rotationally symmetric a-stable process on M". Unlike the Brownian motion case, the explicit 
formula for the density function p{t, x, y) of X with respect to the Lebesgue measure is only known 
for a few special a, such as a = 1. However due to the scaling property of X, one has 

p{t,x,y) = t-"/"p(l,t~^/"x,r^/"y) = t-"/°/(t~i/"(x - y)) for t > and x,y £ M", 

where f{z) is the density function of the symmetric a-stable random variable Xi — Xq in M"'. 
Using Fourier transform, it is not difficult to show (see [9l Theorem 2.1]) that f{z) is a continuous 
strictly positive function on depending on z only through \z\ and that f{z) x Izl"""""^ at infinity. 
Consequently 

1 A , ^ " , 1 on M+ X X M". (1.5) 

\x-y\ 
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In this paper, for two non-negative functions / and g, the notation f ^ g means that there are 
positive constants ci and C2 so that cig{x) < f{x) < C2g{x) in the common domain of definition for / 
and g. For a, 6 G M, af\h := min{a, 6} and aV6 := maxja, h}. However such kind of simple argument 
for (jl.Sp breaks down for the symmetric a-stable-Hke processes on M" when (ajj (x))i<jj<„ = and 
J{x,y) = c{x,y)\x — y\~^~'^ for some a G (0,2) and a symmetric function c{x,y) that is bounded 
between two positive constants in (jl.3p . as in this case, X is no longer a Levy process. 

Two-sided heat kernel estimates for jump processes in M" have only been studied recently. In 
[29], Kolokoltsov obtained two-sided heat kernel estimates for certain stable-like processes in M", 
whose infinitesimal generators are a class of pseudo-differential operators having smooth symbols. 
Bass and Levin [5] used a completely different approach to obtain similar estimates for discrete time 
Markov chain on where the conductance between x and y is comparable to |a; — y|~"~" for a € 
(0,2). In [18j, two-sided heat kernel estimates and a scale-invariant parabolic Harnack inequality 
(PHI in abbreviation) for symmetric a-stable-like processes on d-sets are obtained. Recently in |19j . 
PHI and two-sided heat kernel estimates are even established for non-local operators of variable 
order. Finite range stable-like processes on M"' are studied in |14J. This class of processes is 
very natural in applications where jumps only up to a certain size are allowed. The heat kernel 
estimates obtained in |14] shows finite range stable-like processes behave like discontinuous stable- 
like processes in small scale and behave like Brownian motion in large scale. Processes having 
such properties may be useful in applications. For example, in mathematical finance, it has been 
observed that even though discontinuous stable processes provide better representations of financial 
data than Gaussian processes (cf. pS]), financial data tend to become more Gaussian over a longer 
time-scale (see |31] and the references therein). Our heat kernel estimates in |14j show that finite 
range stable-like processes have this type of property. Moreover, finite range stable- like processes 
avoid large sizes of jumps which can be considered as impossibly huge changes of financial data 
in short time. See [T] for some results on parabolic Harnack inequality and heat kernel estimate 
for more general non-local operators of variable order on M", whose jumping kernel is supported 
on jump size less than or equal to 1. The DeGiorgi-Nash-Moser-Aronson type theory is studied 
very recently in [20j for diffusions with jumps whose infinitesimal generator is of type (jl.3p with 
uniformly elliptic and bounded diffusion matrix (ajj (a;))i<jj<„ and non-degenerate measurable 
jumping kernel J. 

Quite often we need to consider part process of X killed upon exit an open set D C M". 
When X is a Brownian motion, the infinitesimal generator of X^ is the Dirichlet Laplacian ^A^j. 
When X is a rotationally symmetric a-stable process in M", the infinitesimal generator of X^ is 
a Dirichlet fractional Laplacian c A"/-^|£) that satisfies zero exterior condition on D'^. Though the 
transition density function of Brownian motion has been known for quite a long time, due to the 
complication near the boundary, a complete sharp two-sided estimates on the transition density 
of killed Brownian motion in bounded C^'^ domains D (equivalently, the Dirichlet heat kernel) 
have been established only recently in 2002, see [37] and the references therein. Very recently 
in [15], we have obtained sharp two-sided heat kernel estimates for Dirichlet fractional Laplacian 
operator in C^'^ open sets, while in [16] and [T7], we derived sharp two-sided estimates for transition 
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density functions of censored stable processes and relativistic a-stable processes in C ' open sets, 
respectively. 

The rest of the paper is organized as follows. Heat kernel estimates, PHI and a priori Holder 
estimates for stable-like processes and mixed stable-like processes on n-sets in M" is discussed in 
Sections 2 and 3, respectively. In Section 4, we deal with finite range stable- like processes on R", 
while results for diffusions with jumps are surveyed in Section 5. Sections 6 and 7 are devoted to 
sharp heat kernel estimates for symmetric stable processes and censored stable processes in C^'^- 
open sets. To give a glimpse of our approach to the DeGiorgi-Nash-Moser-Aronson type theory 
for non-local operators using probabilistic means, we give an outline of the main ideas in our 
investigation for the following three classes of processes: symmetric stable-like processes on open 
n-sets in M" in Section 2, diffusions with jumps on M" in Section 5 and symmetric stable processes 
in open subsets of M" in Section 6. This paper surveys some recent research that the author is 
involved. See Bass [3j for a survey for related topics on SDEs with jumps, Harnack inequalities and 
Holder continuity of harmonic functions for non-local operators, and Chen [12] for a survey (prior 
to 2000) on potential theory of symmetric stable processes in open sets. 

Throughout this paper, n > 1 is an integer. We denote by m or dx the n-dimensional Lebesgue 
measure in W^, and C^(M") the space of C^-functions on M" with compact support. For a closed 
subset F of M", Cc{F) denotes the space of continuous functions with compact support in F. For 
a Markov process X on a state space E and a subset K C E, we let ax '■= inf{i > : G K} 
and tk '■= ini{t > : Xf ^ K} to denote the first entering and exiting time of K by X. 

Acknowledgement. The author thanks Panki Kim and Renming Song for comments on an 
earlier version of this paper. 



2 Stable-like processes 

A Borel subset F in M" with n > 1 is said to be an n-set if there exist constants ro > 0, C2 > Ci > 
so that 

Ci r" < m{B{x, r)) < C2 r" for all x e F, < r < tq. (2.1) 

In this section and the next, B(x, r) := {y £ F : \x — y\ < r} and | • | is the Euclidean metric in R". 
Every uniformly Lipschitz domain in is an n-set, so is its Euclidean closure. It is easy to check 
that the classical von Koch snowflake domain in R^ is an open 2-set. An n-set can have very rough 
boundary since every n-set with a subset having zero Lebesgue measure removed is still an n-set. 
For a closed n-set F C R" and < a < 2, define 

T = |nEL2(F,m): 

1 f c(x J y \ 
£{u,v) = - {u{x) - u{y)){v{x) - v{y))- '——m{dx)m{dy) (2.3) 

^ Jfxf f ~ y\ 

for u, V £ J^, where c(x, y) is a symmetric function on F x F that is bounded between two strictly 



{u{x) - u{y)y 



FxF 



I n+a 



m{dx)m{dy) < 00 



(2.2) 
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positive constants C4 > C3 > 0, that is, 

C3 < c{x,y) < C4 for m-a.e. x,y £ F. (2-4) 

It is easy to check that {£,J') is a regular Dirichlet form on L^{F,m) and therefore there is 
an associated m-symmetric Hunt process X on F starting from every point in F except for an 
exceptional set that has zero capacity. We call such kind of process a a-stable-like process on F. 
Note that when F = R" and c{x, y) is a constant function, then X is nothing but a rotationally 
symmetric a-stable process on M". 

Theorem 2.1 ([18, Theorem 1.1]) Suppose that F C M" is a closed n-set and < a < 2. Then 
X has a Holder continuous transition density function p{t, x, y) with respect to m. This in particular 
implies that X can he modified to start from every point in F as a Feller process. Moreover, there 
are constants C2 > ci > that depend only on n, a., and the constants Ck, k = 1, - • • ,A in (|2.1|) 
and (j2.4p . respectively, such that 

cimin/t-"/", r-^-— 1 <pit,x,y) < c^mmit'-/'^, —^-—], (2.5) 



for all x,y £ F and <t < 1. 

If F is a global n-set in the sense that (|2.ip holds for every r > 0, then the heat kernel estimates 
in (j2.5p holds for every t > 0. 

Note that in [l8t Theorem 1.1], the dependence of ci,C2 on (Ci,--- ,6*4) in Theorem 12.11 is 
stated for every a except for the case of < a = n < 2. The reason is that in [18j, the on-diagonal 
estimate (Nash's inequality) for the case of a = n < 2 was established by using an interpolation 
method. This restriction can be removed by an alternative way to establish Nash's inequality, see 
[191 Theorem 3.1]. 

The detailed heat kernel estimates such as those in (|2.5p are very useful in the study of sample 
path properties of the processes. For example, the following is proved in jl8j . 



Theorem 2.2 ( |18|, Theorem 1.2]) Under the assumption of Theorem \2.1\ for every x € F, 
F^-a.s., the Hausdorff dimension of X[0, 1] := {Xt : < t < 1} is a An. 

In fact, a much stronger result can be derived from Theorem 12.11 The following uniform 
Hausdorff dimensional result and boundary trace result are established in Remarks 3.10 and 4.4 of 
[6], respectively. 

Theorem 2.3 Let D be an open n-set in M" with n > 2 and X be an a-stable-like process on D. 
Then for every x £ D, 

P^. (dim// X{E) = a dim// E for all Borel sets £; C M+) = 1 

and Tx-a.s. 

dim//(X[0,oo)n9£') =max|l- 

a 
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Here for a time set E C M+, X{E) := {Xt : t G E} and dim//(j4) is the HausdorjJ dimension of a 
set A. 

The approach to Theorem 12. II in [18j is probabiHstic in nature and is motivated by the work of 
Bass and Levin [11[5] on stable-hke processes on Z" and on M". However there are new challenges 
for stable-like processes on n-sets, as paper |3] deals with (possibly non-symmetric) semimartingale 
stable-like processes on R", when restricted to the symmetric processes case, requiring c{x,y) = 
f{x,y — x) and f{x, h) be an even function in /i, while paper [5J is concerned about the transition 
density function estimates for discrete time stable-like Markov chains on Z". 

By Nash's inequality and [H Theorems 3.1 and 3.2], there is a properly exceptional set N <Z F 
and a positive symmetric function p{t, x, y) defined on (0, oo) x [F \ J\f) x (F x J\f) so that x^ y) 
is the density function for Xt under for every x G F \ J\f, 

p{t + s,x,y)= / p{s,x, z)p{t, z,y)m{dz) for every x, y G F \ and t > 
Jf 

and 

p{t, X, y) < c±~'^l'^ for every t > and x,y G F\M . 

Moreover, there is an i5-nest {Fk,k > 1} of compact sets so that M = E \ IJA;>i-^fc ^-iid that for 
every t > and y G F\M, x ^ p{t, x, y) is continuous on each F^. The proof of Theorem 12.11 given 
in [18] relies on the following three key propositions. The first proposition is a tightness result for 
X. 

Proposition 2.4 ([18l Proposition 4.1]) For each ro > 0, ^ > and < B < I, there exists 
< 7 < 1 such that for every < r < ro, 

IPx {tb{x, At) <ir'^) <B for every xeF\J\f. 

Moreover, the constant 7 can be chosen to depend only on {rQ, A, B,n,a) and the constants 
{Ci, 02,03,04) in (gH) and dMD respectively. 

Proposition 2.5 ([18|, Proposition 4.2]) (i) For each a > 0, there exists ci > such that for 
every x £ F\ M, 

(\ d+a 
for every r E (0, 2^-]. (2.6) 

Moreover, the constant ci above can be chosen to depend only on (a, n, a) and on the constants 
{01,02,03,04) in (f2T]l and (i2l4]l . respectively. 

(ii) For each a,b > 0, there exists C2 > such that 

(\ d+a 

for every r E (0, 2^/°^] and such that \x — y\ > br. Moreover, the constant C2 above can be 
chosen to depend only on {a,b,n,a) and on the constants (Ci, C2, C3, C4) in (j2.ip and ()2.4p . 
respectively. 
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The last key proposition is a parabolic Harnack inequality. For this we need to introduce 
space-time process Zg := (V^jX^), where = Vq + s. The filtration generated by Z satisfying 
the usual condition will be denoted as {Ts\ s > 0}. The law of the space-time process s i— > 
starting from (t, x) will be denoted as p(*'^). We say that a non-negative Borel measurable function 
(/(t, x) on [0, co) X F is parabolic in a relatively open subset D of (0, cxd) x F if for every relatively 
compact open subset Di of D, q{t,x) = E^*'"^) {/(Z^-^J for every (t, x) G Di n (0, oo) x (F \ M), 
where tdi = inf{s > : Zg ^ Di}. It is easy to see that for each to > and xq & F \ AA, 
q{t,x) := p{to — t,x,XQ) is parabolic on [0,to) x F- 

For each Rq > 0, we denote ^fi^ := 'y{Ro,l/2,l/2) < 1 the constant in Proposition [23] corre- 
sponding to ro = -Ro and A = B = 1/2. For t < 1 and r < Rq, we define 

QB^{t,x,r) := [t,t + 7/jor"] x {B{x,r) n F \M). 

Proposition 2.6 ([18l Proposition 4.3]) For every Rq > 0, < 6 < ^r^^, there exists c > 
such that for every z£F,0<R<Rq and every non-negative function q on [0, oo) x F that is 
parabolic and bounded on [0,37/jQi?"] x B{z,R), 

sup q{t,y) < c inf q{0,y). 

{t,y)&QRf,{SR°',z,R/3) yeB{z,R/3) 

Moreover, the constant c above can be chosen to depend only on (RQ,6,n,a) and on the constants 
{Ci, €2,03,04) in h2.1\) and \2.4\) respectively. 

Note that the parabolic Harnack inequality implies the elliptic Harnack inequality. 

With the above three propositions, the heat kernel estimates for p(t,x,y) in Theorem 12.11 can 
be established for every t > and x,y £ F\ J\f. That p(t, x, y) is jointly continuous and hence the 
heat kernel estimates hold for every t > and x,y € F comes from the following theorem. 

Theorem 2.7 ( |18l Theorem 4.14]) For every Rq > 0, there is a constant c = c{Rq) > such 
that for every < R < Rq and every bounded parabolic function q in Qjig{0,XQ,max{4:, 4}^'^}R), 

\qis,x) - q{t,y)\ < c\\q\\oo,R R''^ (\t 



|l/a 

holds for {s,x), {t,y) e Qr^{0,Xo,R), where \\q\\oo,R ■= sup(i y)g[o_^^^ max{4,4«}i?°]x(FW) y)l- 
In particular, for the transition density function p{t, x, y) of X, there are constants c > and 
P > such that for any < to < 1, t, s e [to, 2] and {xi,yi) e {F\M) x {F\M) with i = 1,2, 

\p{s,xi,yi) -p{t,X2,y2)\ < Oq^'^^^^^'^ {\t - + \xi - X2I + \yi - 2/2!) • 

Moreover, the constant c above can be chosen to depend only on {Ro,tQ,n,a) and on the constants 
{Oi,C2,03,Oa) in dm) and 
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3 Mixed stable-like processes 



In applications, the stochastic model may have more than one noises. So it is natural to consider 
mixed stable-like processes and a mixture of diffusion and jump-type processes. 

Let -F be a closed global n-set in M". Let (p = be a strictly increasing continuous functions 
on M+, where if) is non-decreasing function on [0, cxd) with ^(r) = 1 for < r < 1 that is either 
the constant function 1 on or there are constants cq > 0, C2 > ci > and 72 > 71 > so that 

cie"'^^ < 4>{t) < C2e^^^ for every 1 < r < 00, (3-1) 

with 

ip{r -|- 1) < co^l'{r) for every r > 1, (3-2) 

and (pi is a strictly increasing function on [0, 00) with (j)i{0) = 0, (/'i(l) = 1 and satisfies the 
following: there exist constants C2 > ci > 0, C3 > 0, and f32 > Pi > such that 

MR) /R\f^2 ^ , ^ 

ci - < < C2 - for every < r < /2 < 00, 3.3) 

\r J vxyT) \r J 
2 

S T 

ds < C3 for every r > 0. (3-4) 



Ms) Mr) 



Remark 3.1 Note that condition (13. 3p is equivalent to the existence of constants 04,05 > 1 and 
Lq > 1 such that for every r > 0, 

C44>i{r) < MLqt) < c^Mf)- 

Denote by d the diagonal ol F x F and J be a symmetric measurable function on F x F \ d 
such that for every {x,y) £ F x F \ d, 

<J(x,y)<i — ^ ^. (3.5) 



\x — yl" (p[C2\x — y\) \x — y|" <p{C4\x — y\ 

For u G L'^{F,m), define := G L'^{F;m) : jpy^p{u{x) — u{y))'^ J{x,y)m{dx)m{dy) < oo| 

and 



£{u,u) := / {u{x) — u{y)){v{x) — v{y))J{x,y)m{dx)m{dy) ior u,v £ J^. (3.6) 
JfxF 

For /? > 0, 

£p{u,u) := £{u,u) + P / u{x)'^m{dx). 

Jf 

It is not difficult to show that {£, J^) is a regular Dirichlet form on L?'{F, m) (see Proposition 
2.2 and Remark 4.10(ii)]. So there is a symmetric Hunt process Y associated with it, starting from 
quasi-every point in F. However the next theorem, which is a special case of [191 Theorem 1.2] (cf. 
[141 Remark 4.4(iv)]), says that X can be refined to start from every point in F. Moreover, it has 
a jointly continuous transition density function p{t, x, y) with respect to the Lebesgue measure on 
F. The inverse function of the strictly increasing function 1 1— > (j)(t) is denoted by (p~^{t). 
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Theorem 3.2 (Theorem 1.2 of |19| ) Under the above conditions, there is a conservative Feller 
process Y associated with {£, T) that starts from every point in F. Moreover the process Y has a 
continuous transition density function on (0, oo) x F x F with respect to the measure m, which has 
the following estimates. There are positive constants ci > 0, C2 > and C > 1 such that 



^{t)'^ \x — y\"' (j){ci\x — y\) J ' ' \(/) ^(t)" \x — y\"-(l){c2\x — y\ 

for every t G (0, 1] and x,y £ F. Moreover, when 1^=1, the above heat kernel estimates hold for 
every t > and x,y £ F. 

We now give some examples such that Theorem 13.21 appHes. 

Example 3.3 If there is < ai < 02 < 2 and a probabihty measure on [01,02] such that 

0(r) := ( / \~''u{da) 



\-Jai 

then conditions (|3.3p - (|3.4p are satisfied with ip = 1. Clearly, i;^ is a continuous strictly increasing 
function with i;^(0) = and = 1. The condition (j3.3p is satisfied with 7 = 1 because 

1 (t>(r) 1 
< — — r < for any r > 0. 

2"! - 4>{2r) - 2"2 

For r > 0, by Fubini's theorem, 

l-r i-r 1-02 pa2 i i 2 

/ ih'^' = / / ^'~''^^da)ds = / r2-^(da) < -— 

Jo (P{s) Jo Jqi 7ai 2 - a 2 - Q2 0(r) 

and so condition ()3.4p is satisfied. In this case, 

A particular case is when is a discrete measure. For example, v is a discrete measure concentrate 
on a, (3 G (0,2). In this case, J{x,y) = ^^^''y^nla + jJi^irFF' where Ci{x,y) are two symmetric 
functions that are bounded between two positive constants, and 



mm < r", r^ 



}, 0-i(r) X maxjr^/", r^/^}. 



Theorem 13.21 gives the precise heat kernel estimates for mixed stable- like processes on F. When 
F = M", Theorem 13.21 in particular gives the heat kernel estimate for Levy processes on which 
are linear combinations of independent symmetric o-stable processes. Of course. Theorem 13.21 holds 
much more generally, even in the case of F = M". □ 



Example 3.4 Let Y = {Yt,t > 0} be the relativistic a-stable processes on M" with mass tuq > 0. 
That is, {Yt,t > 0} is a Levy process on R"" with 

E[eMi{tYt - Yo))] = exp (t (m^ - del' + m^)"/^)) . 
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where a E (0,2). It is shown in [22] that the corresponding jumping intensity satisfies 

^(mo|x - y\) 



J{x,y) 



where ^{r) x e~''(l + near r = cxd, and ^'(r) = 1 + ^"(0)r^/2 + o(r^) near r = 0. So 

conditions (j3.ip - ()3.4p ) are satisfied with 71 > for the jumping intensity kernel for every relativistic 
a-stable processes on R". 

When a = 1, the process is called a relativistic Hamiltonian process. In this case, the heat 
kernel can be written as 

p{t,x,y) = , ^ [ e™°*e-V(l-?/P+*^)(l-P+-g)dz, 

see |24] . For simplicity, take niQ = 1. It can be shown that for every t > and (x,y) G M" x R", 

(1 V (|x -y\+ tf^^^ e~^' Vli-.P+t^ 



< p{t,x,y) 



(|x-y|+t)"+i . 

This in particular implies that for every fixed to > 0, there exist ci, • • • , C4 > which depend on to 
such that 

for every t £ (0, to] ^-nd x,y £ M", which is a special case of Theorem 13. 2i □ 



The following construction of Meyer |32J for jump processes played an important role in our 
approach in [19]. Suppose we have two jump intensity kernels J{x,y) and Jo{x,y) on F x F such 
that their corresponding pure jump Dirichlet forms given in terms of (13. 6p with = T>{£) are 
regular on F. Let Y = {Yt,t > 0, P^, x e F \ AA} and = {Y^^\t > 0, P^,x G F\Mq} be 
the processes corresponding to the Dirichlet forms whose Levy densities are J{x,y) and Jo(x,y), 
respectively. Here M and A/q are the properly exceptional sets of Y and Y^^^ , respectively. Suppose 
that jQ{x,y) < J{x,y) and 



:= / (J(x,y) - Joix,y))m{dy) < c, 



F 



for all X G F. Let 



Ji{x,y) := J{x,y) - Jo{x,y) and = "^^^^'f^ (3.7) 

J(x) 

Then we can construct a process 1" corresponding to the jump kernel J from 1"^'^) as follows. Let 
5i be an exponential random variable of parameter 1 independent of Y^^\ let Ct = Jq J(Yj^^) ds, 
and let Ui he the first time that Ct exceeds Si. We let 1^ = Yj^^ for < s < Ui. 
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At time Ui we introduce a jump from Yiy^- to Zi, where Zi is chosen at random according to 
the distribution q{Yi/-^-, y). We set Yui = Zi, and repeat, using an independent exponential S2, 
etc. Since J^{x) is bounded, only finitely many new jumps are introduced in any bounded time 
interval. In [32] it is proved that the resulting process corresponds to the kernel J. See also [26] . 
Note that if A/q is the properly exceptional set corresponding to Y^^^ , then this construction gives 
that the properly exceptional set M for Y can be chosen to be a subset of A/q. 

Conversely, we can also remove a finite number of jumps from a process Y to obtain a new 
process Y^^\ For simplicity, assume that Jo{x,y)Ji{x,y) = 0. Suppose one starts with the process 
Y (associated with J), runs it until the stopping time Si = mi{t : Ji{Yt-,Yt) > 0}, and at that 
time restarts Y at the point Is^-. Suppose one then repeats this procedure over and over. Meyer 
|32j proves that the resulting process Y^'^^ will correspond to the jump kernel Jq. In this case 

Assume that the processes Y and Y^^^ have transition density functions p{t, x, y) and p^'^^ (t, x, y), 
respectively. Let {Tt}t>Q be the filtration generated by the process Y^^\ The following lemma is 
shown in [U Lemma 2.4] and in ^ Lemma 3.2]. 

Lemma 3.5 (i) For any A ^ J^t, 

({Ys = yP for allQ<s<t]r\A^ > e~*ll^ll°° P^(^). 

(ii) //||Ji||oo < 00, then 

p{t,x,y) < p'-^\t,x,y) +t||Ji||oo- 

The use of Lemma 13.51 can also significantly simplify the proofs in [TB] for results in the last 
section. The relation between Y and Yq can be viewed as the probabilistic counterpart of the 
Trotter's semigroup perturbation method. For example, the proof for Propositions 12.5112.61 can be 
simplified by using Lemma l3.5i See the proof of Propositions 4.9 and 4.11 of [19j in this regard. 

Comparing with Theorem 12. H Theorem 13.21 says that the rate function for stable processes 
of mixed type associated with (|3.5p - (j3.6p is cp. Parabolic Harnack inequality and a prior Holder 
estimate also hold for parabolic functions of X, with this rate function (j). For each r,t > 0, we 
define 

Q{t, X, r) := [t, t + 70(r)] x {B{x, r) D F) . 

Theorem 3.6 ([19', Theorem 4.12]) For every < (5 < 7, there exists ci > such that for every 
z £ F , R & (0, 1] (resp. R > when 71 = 72 = and every non-negative function h on [0, cxd) x F 
that is parabolic and bounded on [O,70(2i?)] x B{z,2R), 

sup h{t,y)<ci inf h{0,y). 

In particular, the following holds for t < 1 (resp. t > when 71 = 72 = Oj. 

sup p{s,x,y)<c inf p{{l + j)t,x,y). (3.8) 

{s,y)GQ({l-7)i,2,0-i{t)) y(iB(z,<f>~^(t)) 
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Proposition 3.7 ([19], Proposition 4.14]) For every Rq G (0,1] (resp. Rq > when 71 = 72 = 
0), there are constants c = c{Ro) > and k > such that for every < R < Rq and every bounded 
parabolic function h in Q{0,xq, 2R), 

\h{s,x) - h{t,y)\ < c||/i||oo,Fi?-" {rWt-s\) + p{x,y)Y 

holds for {s,x), {t,y) £ Q{0,xo,R), where ||/i||oo,F := sup(t^^)g[o,^0(2j?)]xF ^"^ particular, 

for the transition density function p{t,x,y) of X, for any to £ (0,1) (resp. any T > and any 
to £ (0, T) when 71 = 72 = Oj, there are constants c = c{to) > and k > such that for any 
t, s G [to, 1] (resp. t,s £ [to, T]) and (xj, yi) £ F x F with i = 1,2, 

\p{s, xi,yi)- p{t, X2,y2)\ < c ('^"Hl* + p{xi, X2) + p{yi,y2)Y ■ 

4 Finite range stable-like processes 

A finite range a-stable-like process X on M" is a symmetric Hunt process on M" of purely discon- 
tinuous type whose jumping kernel is J{x,y) = j^l^^jn+a '^[\x-y\<K.}-, where a £ (0,2), k > and 
c{x, y) is a symmetric function on M" x M" that is bounded between two positive constants. The 
Dirichlet form {£,J-) associated with X on L?'(W^,m) is given by 



r = \u£L\W^;m): I (^(x)^ u{y)f i^^^_^^^^^rn{dx)m{dy) < 00 \ (4.1 



^xR'^ 1-^ ~ y\ 



nGL2(M";m): / M^L_^iMl!m(dx)m(dy) < 00 [> , 



1 f c(x i/) 
£{u,v) = - {u{x) - u{y)){v{x) - v{y))- '-T:;—^l{\x-y\<K}T^{dx)m[dy) (4.2) 

^ J FxF F ~ y\ 

for u,v £ J-. The L^-infinitesimal generator of X and {£,J-) is a non-local (integro-differential) 
operators £ on M" with measurable symmetric kernel J{x,y) = ^J^y[^\a 'i-{\x-y\<K}- 

Cu{x) = lim / {u{y) — u{x))J{x,y)dy. 

ej,0 J|„gK": \y-x\>e'\ 



■ \y-x\>e} 

Theorem 4.1 ( |14L Proposition 2.1 and Theorems 2.3 and 3.6]) The finite range stable-like 
process X has a jointly continuous transition density function p{t, x, y) and so X can be refined to 
start from every point on M". Moreover the following sharp two-sided heat kernel estimates hold. 

(i) There is i?* £ (0, 1) so that for every t £ (0, -R"] and x,y £W^ with \x — y\ < R^ 

(ii) There exists C=k £ (0, 1) such that for x,y £ M" with \x — y\ > max{t/C*, i?*}, 

J. \ c\x—y\ / I I 

^ \ I \ M F ~ y| 



p{t,x,y)^(- r = exp -c|x - 7/1 log 



x-y\J V t 
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(iii) For t > R'^ orx,y £W with \x - y\ £ [R^,t/C^], 

n/2 ( C\x-y? 



The following weighted Poincare inequality for non-local operators together with Lemma [37 
played a crucial role in our proof of Theorem 14.11 in [14] . In the remainder of this paper, B{x,r) 
denotes the Euclidean ball in M" with radius r centered at x. 

Theorem 4.2 ( [14. Proposition 3.2]) Suppose that J{x,y) is a symmetric non-negative kernel 
on M" X M" such that J{x, y) = when |x — y| > 1 and 



Ki\x - y\ 



" " < J{x,y) < K2\x — y\ ^ when \x — y\ < 1 



for some constants ki,K2 > and < a < P < 2. Let (j){x) := c (l — 1b(o,i)(^)) 
where c > is the normalizing constant so that (f){x)dx = 1. Then there is a positive constant 
ci = ci(n, a,/3) independent of r > 1, such that for every u £ L^{B{0,l),4>dx), 

/ {u{x) — u^)'^(p{x)dx < ci / {u{x) — u{y))'^ r"'^'^J{rx, ry)\/ (j){x)(t){y) dxdy. 

Jb{o,i) Jb{0,1)xB{0,1) 



Here 



'^4> ■= Ib{o,i) u{x)(l){x)dx. 



5 Diffusions with jumps 

In this section, we consider symmetric Markov processes on that have both the diffusion and 
pure jumping components. More precisely, consider the following regular Dirichlet form (<?,J^) on 
L^(R"';m) given by 



£{u^v) = - I Vu{x) ■ A{x)Vv{x)dx + I {u{x) — u{y)){v{x) — v{y))J{x,y)dxdy, 



2 



(5.1) 



where A(x) = {aij{x))i<ij<n is a measurable n x n matrix- valued function on M" that is uniform 
elliptic and bounded in the sense that there exists a constant c > 1 such that 

n n n 

4' < J2 (^)^*^^- ^ ^ E ^^^^y ^, (6 , • • • , ^d) e M^ (5.2) 

1=1 *j=i *=i 

and J is a symmetric non-negative measurable kernel on R" x R" such that there are positive 
constants kq > 0, and (3 E (0, 2) so that 

J(x,y) < Ko|x - for|x-y|<5o, (5.3) 

and that 

sup / {\x — yl"^ A 1) J{x,y) dy < oo. (5.4) 
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Clearly under condition (j5.3p . condition (j5.4p is equivalent to 

sup / J{x,y) dy < CO. 

x&R" J{y€M.^^:\y-x\>l} 

By the Dirichlet form theory, there is an M'^-valued symmetric Hunt process X associated with 
{£,J-). The L^-infinitesimal generator of X ia a non-local (pseudo-differential) operators C on M": 

1 " 9 f du{x)\ f 
Cu{x) = o 1^ (^iji^)-^ -hlim / (n(y) - u{x)).J{x,y)dy, (5.5) 

When the jumping kernel J = in (j5.5p and (jS.ip . £ is a uniform elliptic operator of divergence 
form and X is a symmetric diffusion on M". It is well-known that X has a joint Holder continuous 
transition density function p{t,x,y), which enjoys the celebrated Aronson's two-sided heat kernel 
estimate (II. ip . 

When A{x) = in (jS.ip and J is given by 

J{x,y)-^- |^-^T| (5.6) 

\x - y\'^(l)[\x - y\) 

where (j) a strictly increasing continuous function (p : M_j_ — > R_(_ with 0(0) = 0, and = 1 that 
satisfies the conditions p.3p - ()3.4p with (p in place of cpi there, the corresponding process X is a 
mixed stable-like process on appeared in the previous section. We know from Theorem 13.21 that 
there are positive constants < ci < C2 so that 

cip'it, \x - y\) < p{t,x,y) < C2P^{t, \x - y\) for t>0,x,y £ M", 

where 

I^it,r):= ( ^-\t)-- A (5.7) 



r"-(p[r) ^ 

with </)~^ being the inverse function of (p. 

In this section, we consider the case where both A and J are non-trivial in (j5.5p and (jS.ip . 
Clearly such non-local operators and diffusions with jumps take up an important place both in 
theory and in applications. However there are very limited work in literature for this mixture case 
on the topics of this paper until very recently. One of the difficulties in obtaining fine properties 
for such an operator C and process X is that they exhibit different scales: the diffusion part has 
Brownian scaling r i— > while the pure jump part has a different type of scaling. Neverthe- 
less, there is a folklore which says that with the presence of the diffusion part corresponding to 
^ Yll j=i ^ af") ' better results can be expected under weaker assumptions on the jumping 

kernel J as the diffusion part helps to smooth things out. Our investigation in [20] confirms such 
an intuition. In fact we can establish a priori Holder estimate and parabolic Harnack inequality 
under weaker conditions than (15. 6p . We now present the main results of [20]. Let W^^'^{W) denote 
the Sobolev space of order (1,2) on M"; that is, W^^'^{W) := {/ G L2(M";m) : V/ e L2(M";m)}. 
It is not difficult (see Proposition 1.1 of [20]) to show that under the conditions (j5.2p - (j5.4p . the 
domain of the Dirichlet form of (jS.ip is characterized by 
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and that (Theorem 2.2 of [20j) the corresponding process X has infinite hfetime. Let Z = {Zt := 
(Vo — t,Xt),t > 0} denote the space-time process of X. We say that a non-negative real valued 
Borel measurable function h{t, x) on [0, oo) x M" is parabolic (or caloric) on D = (a, b) x B{xQ,r) if 
there is a properly exceptional set M C M" such that for every relatively compact open subset Di 
oiD, 

Mt,x)=E(*'-)[/i(Z,^J] 

for every {t,x) e Di D ([0,cx)) x (M" \Af)), where td^ = inf{s > : ^ Di}. We remark 
that in Sections 2 and 3 the space-time process is defined to be (Vq -|- t, Xt) but this is merely 
a notational difference. (For reader's convenience, we keep the notations same as those in the 
references [M [El 120].) 

Theorem 5.1 (Theorem 1.2 of [20]) Assume that the Dirichlet form given by (jS.ip sat- 

isfies the conditions (|5.2|) - (j5.4|) and that for every < r < 5q, 

inf inf / J(x, z)dz > 0. (5.8) 

x(iB{xo,r/l&) JB{yo,r/16) 

Then for every Rq £ (0,1], there are constants c = c{Ro) > and n > such that for every 
< R < Rq and every bounded parabolic function h in Q{0, Xo,2R) := (0, 4i?^) x B{xo,2R), 

\h{s,x)-h{t,y)\ <c\\h\\^,RR-^ (\t-s\^/^ + \x-y\y 

holds for {s,x), {t,y) E {3R'^,4R'^) x B{xq,R), where \\h\\oo,R ■= s^V{t,y)&[o,AR^]y.M.^\j^\Kt^y)\- In 
particular, X has a jointly continuous transition density function p{t, x, y) with respect to the 
Lebesgue measure. Moreover, for every to G (0, 1) there are constants c > and k > such 
that for any t, s G {tQ, 1] and {xi,yi) e x with i = 1,2, 

\p{s,xi,yi) -p(t,X2,y2)\ < ctg '-"^'^-'^^ (^\t - s|^/^ + \xi - X2I + \yi - 2/2I) ■ 

In addition to (j5.2|) - (|5.4p and (j5.8p . if there is a constant c > such that 

J{x,y) < — J{z,y)dz whenever r < ||x — y| A 1, a;, y G M", (5.9) 

^" Jb(x,t) 

then the following parabolic Harnack principle holds for non-negative parabolic functions of X. 

Theorem 5.2 (Theorem 1.3 of [20| ) Suppose that the Dirichlet form {£,!F) given by (j5.ip sat- 
isfies the condition (|5.2p - (|5.4p . (|5.8p and (|5.9p . For every 5 € (0, 1), there exist constants c\ = ci{5) 
and C2 = C2{5) > such that for every z £ M", to ^ 0, < i? < ci and every non-negative function 
u on [0,00) X M" that is parabolic on {tQ,to + 66R'^) x B{z,AR), 

sup u{ti,yi)<C2 inf u{t2,y2), (5.10) 
(ti,yi)eQ- (t2,y2)eQ+ 

where Q- = (to + to + 25R^) x B{xq, R) and Q+ = (to + 35i?2, to + 'i5R'^) x B{xo, R). 
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We next present a two-sided heat kernel estimate for X when J(x, y) satisfies the condition 
(|5^ . Clearly (fO]) - (fO) . (f5^ and (fOI) are satisfied when (f5^ holds. Recall that functions 
p'^{t,x,y) and {t, x,y) are defined by (|1.2p and (|5.7p . respectively. 

Theorem 5.3 (Theorem 1.4 of [20j) Suppose that (j5.2p /loWs and i/iai i/ie jumping kernel J 
of the Dirichlet form given by (jS.ip satisfies the condition (|5.6p . Denote by p{t,x,y) the 

continuous transition density function of the symmetric Hunt process X associated with the regular 
Dirichlet form of (jS.ip with the jumping kernel J given by (|5.6p . There are positive constants 

Ci, i = 1, 2, 3, 4 such that for every t > and x, y E M", 

ci A <p-\t)-'') A (/(t,C2|x - y\)+pi{t, \x - y\)) 

< pit,x,y)<C3 (t'^/^ArHt)-'^)/\{p''it,c4x-y\)+pi{t,\x-y\)). (5.11) 

Here p'^ and p^ are the functions given by (jl.2p and ()5.7p . respectively. 

When A(x) = Inxm the nxn identity matrix, and J(x, y) = c\x — y|~"^" for some a E (0, 2) in 
(jS.ip . that is, when X is the independent sum of a Brownian motion W on M" and an isotropically 
symmetric a-stable process Y on M", the transition density function p{t,x,y) can be expressed 
as the convolution of the transition density functions of W and Y, whose two-sided estimates 
are known. In [35], heat kernel estimates for this Levy process X are carried out by computing 
the convolution and the estimates are given in a form that depends on which region the point 
{t, X, y) falls into. Subsequently, the parabolic Harnack inequality (j5.10p for such a Levy process 
X is derived in |35j by using the two-sided heat kernel estimate. Clearly such an approach is 
not applicable in our setting even when </>(r) = r°', since in our case, the diffusion and jumping 
part of X are typically not independent. The two-sided estimate in this simple form of (|5.1ip is 
a new observation of [20] even in the independent sum of a Brownian motion and an isotropically 
symmetric a-stable process case considered in [35]. 

The approach in [20j employs methods from both probability theory and analysis, but it is 
mainly probabilistic. It uses some ideas previously developed in [Ul^l dSl dHl dl]- To get a priori 
Holder estimates for parabolic functions of X, we establish the following three key ingredients. 

(i) Exit time upper bound estimate: 

^x[TB{xo,r)] < cir^ for X E B{xo,r), 
where TB{xo,r) ■= inf{t > : Xt ^ B{xQ,r)} is the first exit time from B{xo,r) by X. 

(ii) Hitting probability estimate: 

2 

^ B{x, s)) < for every r E (0, 1] and s > 2r. 
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(iii) Hitting probability estimate for space-time process Zt = (Vq — t^Xf): for every x G M", 
r G (0, 1] and any compact subset A C Q{x,r) := (0, r^) x B{x,r), 

'{cTA<rr)>C3 ^^^^ , 

where by slightly abusing the notation, aA '■= {t > : Zf G A} is the first hitting time of 
A, Tr := inf{t > : Zf ^ Q{x,r)} is the first exit time from Q{x,r) by Z and run+i is the 
Lebesgue measure on M"^"-^. 

Here we use the following notations. The probability law of the process X starting from x is 
denoted as and the mathematical expectation under it is denoted as E^;, while probability law 
of the space-time process Z = {V,X) starting from {t,x), i.e. {Vo,Xo) = {t,x), is denoted as p(*'^) 
and the mathematical expectation under it is denoted as E^*'^). To establish parabolic Harnack 
inequality, we need in addition the following. 

(iv) Short time near-diagonal heat kernel estimate: for every to > 0, there is C4 = C4(to) > such 
that for every xq G M" and t £ (0, to]; 

> C4t-"/2 for x,y £ B{xo,Vt/2). 

Here p^(^o,Vi) ^^j-^g transition density function for the part process x^^^^'^'^ of X killed 
upon leaving the ball B{xo, \/t). 

(v) Let i? < 1 and S < 1. Qi = [to + 26R'^/3, to + 55R^] x B{xo, 3R/2), Q2 = [to + 5R^/3, to + 

115R'^/2] X B{xo,2R) and define Q_ and Q+ as in Theorem Let h : [0,oo) x M" ^ R+ 
be bounded and supported in [0, 00) x B{xo, 3Ry. Then there exists C5 = 05(6) > such that 

E(*1'^i)[/i(Z,qJ] < C5E(*2'^2)[/i(Z,qJ] for G Q_ and (t2,y2) G Q+. 

The proof of (iv) uses ideas from p], where a similar inequality is established for finite range 
pure jump process. However, some difficulties arise due to the presence of the diffusion part. 

The upper bound heat kernel estimate in Theorem 15. 31 is established by using method of scaling, 
by Meyer's construction of the process X based on finite range process X^^'), where the jumping 
kernel J is replaced by J(x, y)t{\x-y\<x}, and by Davies' method from [Hj to derive an upper bound 
estimate for the transition density function of X^^^ through carefully chosen testing functions. Here 
we need to select the value of A in a very careful way that depends on the values of t and [x — y[. 

To get the lower bound heat kernel estimate in Theorem \5.3\ we need a full scale parabolic 
Harnack principle that extends Theorem 15.21 to all ii > with the scale function (/>(i?) := R'^ A(j){R) 
in place oi R ^ R? there. To establish such a full scale parabolic Harnack principle, we show the 
following. 

(iii') Strengthened version of (iii): for every x G M", r > and any compact subset A C 
Q(0,x,r) := [O,7o0(r)] xS(x,r), 
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Here 70 denotes the constant 7(1/2, 1/2) in Proposition 6.2 of [20] . 

(vi) For every 5 £ (0,70], there is a constant cq = 06(7) so that for every < R < 1, r £ (0, i2/4] 
and {t, x) G Q(0, z, R/3) with < t < 7o^(-R/3) - 54>{r), 



C6- 



where U{t,x,r) := {t} x B{x,r). 

With the full scale parabolic Harnack inequality, the lower bound heat kernel estimate can then be 
derived once the following estimate is obtained. 

(vii) Tightness result: there are constants cy > 2 and > such that for every t > and 
X, y G M" with \x — y\ > CTCpit), 

X,eB{y,crrHt)))>cs '^^''^''^^^ 



X - yY'(l)(\x - y\) 

6 Dirichlet heat kernel estimates for symmetric stable processes 

Many times one encounters part process of X killed upon exiting a open set D. The infinitesimal 
generator of X^ is the infinitesimal generator £ of X satisfying Dirichlet boundary or zero 
exterior condition. It is a fundamental problem both in analysis and in probability theory to 
study precise estimate for the transition density function of X^ (or equivalently, the Dirichlet heat 
kernel of C^). However due to the complication near the boundary, two-sided estimates on the 
transition density of killed Brownian motion in bounded C^'^ domains D (equivalently, the Dirichlet 
heat kernel) have been established only recently in 2002; see |37) and the references therein. In this 
section, we survey the recent result from |15) on sharp two-sided estimates on the transition density 
function p£i{t,x,y) of part process X^ of a rotationally symmetric a-stable process killed upon 
leaving a C^'^ open set D. The infinitesimal generator of X^ is the fractional Laplacian cA"/^|£) 
satisfying zero exterior condition on D'^. 

Recall that an open set D in M" (when n > 2) is said to be a C^'^ open set if there exist a 
localization radius i?o > and a constant Aq > such that for every z £ dD, there is a C^'^-function 
= 02 : M"-i M satisfying (/>(0) = V</>(0) = 0, ||V(/>||oo < Aq, \V(I){x) - V(/>(z)| < Ao|x - z\, and 
an orthonormal coordinate system CSz- y = (yi, • • • ,yn-i,yn) '■= (y, Vn) with its origin at z such 
that 

B{z, Ro)nD = {yG B{0, Rq) : y„ > 0(y)} , 

where the ball B{0,Rq) on the right hand side is in the coordinate system CSz- The pair {Rq,Aq) 
is called the characteristics of the C^'^ open set D. We remark that in some literatures, the C^'^ 
open set defined above is called a uniform C^'^ open set as {Rq,Aq) is universal for every z £ dD. 
For X S M", let 6d{x) denote the Euclidean distance between x and W^. By a C^'^ open set in M 
we mean an open set which can be written as the union of disjoint intervals so that the minimum 
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of the lengths of all these intervals is positive and the minimum of the distances between these 
intervals is positive. Note that a C^'^ open set can be unbounded and disconnected. 

Theorem 6.1 ( [15|, Theorem 1.1]) Let D be a C^'^ open subset of MP with n > 1 and Sd{x) the 
Euclidean distance between x and D'^. 



(i) For every T > 0, on (0, T]x D x D, 

n+a 

I I t^l"" \ 

/ . \ . — n /rv / - . *' 1 

PD 



(ii) Suppose in addition that D is bounded. For every T > 0, there are positive constants c\ < C2 
so that on [T, oo) x D x D, 

cie-'-^'doixT'^SDiyT'^ < PD{t,x,y) < C2e-'^^HD{xr'^ doivT'^ 
where Ai > is the smallest eigenvalue of the Dirichlet fractional Laplacian (— A)"/^|£). 

By integrating the two-sided heat kernel estimates in Theorem 16.11 with respect to t, one can 
easily recover the following estimate of the Green function G£){x,y) = pD{t, x,y)dt, initially 
obtained independently in [21] and [30] when n > 2. 



Corollary 6.2 ( [15. Corollary 1.2]) Let D be a bounded C^'^-open set in M" with n> 1. Then 
on D X D, 



GD{x,y) 



F^(lA ^-^-;y^ ) whenn>a, 
[(5z.(x)fe(y))("-^)/^AM£):gjM^ whenr^ 



I < a. 



Theorem 16. iT i) is established in [15] through Theorems 16. 31 and 16. 41 which give the upper bound 
and lower bound estimates, respectively. Theorem I6.1l fii) is an easy consequence of the intrinsic 
ultracontractivity of the symmetric a-stable process in a bounded C^'^ open set. In fact, the upper 
bound estimates in both Theorem 16.11 and Corollary 16.21 hold for any domain D with (a weak 
version of) the uniform exterior ball condition in place of the C^'^ condition, while the lower bound 
estimates in both Theorem 16. II and Corollarv 16.21 hold for any domain D with the uniform interior 
ball condition in place of the C^'^ condition. 

We say that D is an open set satisfying (a weak version of) the uniform exterior ball condition 
with radius tq > if for every z S dD and r G (0,ro), there is a ball of radius r such that 
B"" CW\D and dB^ ndD = {z}. 

Theorem 6.3 ( [151, Theorem 2.4]) Let D be an open set inM" that satisfies the uniform exterior 
ball condition with radius ro > 0. For every T > 0, there exists a positive constant c = c{T,rQ,a) 
such that for t £ (0, T] and x,y £ D, 

PD{t,x,y)<ct-^l''\l^-—^ lA^ . (6.1) 
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An open set D is said to satisfy the uniform interior ball condition with radius ro > in the 
following sense: For every x £ D with 6d{x) < ro, there is Zx G dD so that \x — z^l = Snix) and 
B{xo,rQ) C D for xq := Zx + ^0(2; — Zx)/\x — Zx\ ■ It is well-known that any (uniform) C^'^ open 
set D satisfies both the uniform interior ball condition and the uniform exterior ball condition. 

Theorem 6.4 ( |15|. Theorem 3.1]) Assume that D is an open set in satisfying the uniform 
interior ball condition. Then for every T > there exists a positive constant c = c{rQ,a,T) such 
that for all {t, x, y) € (0, T] x L> x D, 

There are significant differences between obtaining two-sided Dirichlet heat kernel estimates for 
the Laplacian and the fractional Laplacian, as the latter is a non-local operator. Our approach in 
|15j is mainly probabilistic. It uses only the following five ingredients: 

(i) the upper bound heat kernel estimate for the rotationally symmetric a-stable process X in M" 

and the stable-scaling property of X; 

(ii) the Levy system of X that describes how the process jumps; 

(iii) the mean exit time estimates from annuli and from balls; 

(iv) the boundary Harnack inequality of X in annuli (when n>2) and in intervals (when n = 1), 
and the parabolic Harnack inequality of X; 

(v) the intrinsic ultracontractivity of X in bounded open sets. 

The upper bound heat kernel estimate of X on gives an upper bound for pD{t,x,y), while 
the Levy system is the basic tool used throughout our argument as the symmetric stable process 
moves by "pure jumping". To get the boundary decay rate of pDit,x,y), we use the boundary 
Harnack inequality and the domain monotonicity of the killed stable process X^ in D by comparing 
it with certain truncated exterior balls (i.e. annulus) as well as interior balls. The mean exit 
time estimate for an annulus is applied with the help of the boundary Harnack inequality to 
get the boundary decay rate in the upper bound heat kernel estimates. The two-sided estimates 
in the ball B = B{0,1): E^[rB] x 6b{x)'^^'^ is used to get the two-sided estimate on the first 
eigenfunction in balls. The latter is then used to get the boundary decay rate for the lower bound 
estimate in pD{t, x, y). The parabolic Harnack inequality allows us to get pointwise lower bound on 
PD{t,x,y) from the integral of 1— > pD{t/2,x,w) over some suitable region. When X^ is intrinsic 
ultracontractive, pD{t,x,y) is comparable to ct4>D{x)4>D{y) for some ct > and a good control 
is known for q when t is above a certain large t^, where </>_d is the positive first eigenfunction of 
(— A)°/^|£), the infinitesimal generator of X^ . 
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Note that the large time heat kernel estimate in Theorem 16. If ii) requires D to be bounded. See 
|23j for recent results on large time sharp heat kernel estimates for symmetric stable processes in 
certain unbounded C^'^ open sets. 

The approach developed in [15] is quite general in principle and can be adapted to study heat 
kernel estimates for other types of jump processes in open subsets and their perturbations, such as 
censored stable processes to be discussed in next section. 

7 Dirichlet heat kernel estimates for censored stable processes 

Censored a-stable processes in an open subset of M" were introduced and studied by Bogdan, 
Burdzy and Chen in [8j. Fix an open set D in M" with n > 1. Define a bilinear form £ on C^{D) 
by 

£{u,v):=\f f {u{x)-u{y)){v{x)-v{y))- ''—^dxdy, u,veC^{D), (7.1) 

^ Jd Jd f ~ y\ 

where c > is a constant. Using Fatou's lemma, it is easy to check that the bilinear form {£, C^{D)) 
is closable in L'^{D,dx). Let J- be the closure of C^{D) under the Hilbert inner product £i := 
+ ( ■ ) ■ )L'^{D,dx)- is noted in [8j, {£,J^) is Markovian and hence a regular symmetric Dirichlet 
form on L'^{D,dx), and therefore there is an associated symmetric Hunt process X = {Xt,t > 
0,Fx,x G D} taking values in D. The process X is called a censored (or resurrected) a-stable 
process in D. 

Let y be a rotationally symmetric a-stable process in M" with jumping kernel c\x — yl"""". 
For any open subset D of M", we use to denote the subprocess of Y killed upon exiting from 
D. The following result gives two other ways of constructing a censored a-stable process. 

Theorem 7.1 ((S], Theorem 2.1 and Remark 2.4]) The following processes have the same dis- 
tribution: 

(i) the symmetric Hunt process X associated with the regular symmetric Dirichlet form {£,J^) on 

L'^{D,dx); 

(ii) the strong Markov process X obtained from the killed symmetric a-stable-like process Y^ in 
D through the Ikeda-Nagasawa-Watanabe piecing together procedure; 

(Hi) the process X obtained from Y^ through the Feynman-Kac transform e-^o i^D(yP)ds y^j^^f^ 

■■= / ] '-r^ni^dy. 

J DC \x - 

The Ikeda-Nagasawa-Watanabe piecing together procedure mentioned in (ii) goes as follows. 
Let Xt{Lo) = Y/^{lj) for t < td{uj). If Y^^_{lo) ^ D, set Xt{Lo) = d for t > tb(c^). If Y^^^{lu) G D, 
let XT-jy{u!) = Y^_{uj) and glue an independent copy of Y^ starting from YJ^_{uj) to Xrjj{uj). 
Iterating this procedure countably many times, we obtain a process on D which is a version of the 
strong Markov process X; the procedure works for every starting point in D. 
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For any open n-set D in MJ^, define 



and 



'D Jd f — yi 



f 1 f f C f 

{u,v) := - {u{x) - u{y)){v{x) - v{y))-r——r^dxdy, u,v e T''" . 



As we see from Section 2, the bilinear form (£^^^ is a regular symmetric Dirichlet form on 

L'^{D,dx). The process X on D associated with (if^^^, JF"^*^^) is called in [8j a reflected a-stable 
process on D. By Theorem 12. 11 X has a Holder continuous transition density function p{t, x, y) on 
(0, oo) X D X D and for every Tq > 0, there are positive constants ci, C2 so that for t £ (0, Tq] and 

lA^^l <p(t,x,y) <C2t-"/" ( lA^^ 1 . (7.2) 

\ \x-y\J 

This in particular implies that X can start from every point in D. When D is an open n-set in M", 
the censored a-stable-like process X can be realized as a subprocess of X killed upon leaving D 
(see [SI Remark 2.1]). It is proved in [5] that when L> is a global Lipschitz domain and a £ (0, 1], 
then X and X are the same and so X has a sharp two-sided heat kernel estimate (|7.2p in this case. 
Hence in the following we will concentrate on the case of a G (1, 2). The next theorem gives a sharp 
two-sided heat kernel estimate for the transition density function p£){t,x,y) of censored a-stable 
process in an C^'^ open set with a £ (1, 2). 

Theorem 7.2 ( |16L Theorem 1.1]) Suppose that n > 1, a £ (1,2) and D is a C^'^ open subset 
o/M"". Let S^ix) be the Euclidean distance between x and D^. 

(i) For every T > 0, on {0,T] x D x D 

(ii) Suppose in addition that D is bounded. For every T > 0, there exist positive constants ci < C2 
such that for all {t,x,y) £ [T, oo) x D x D, 

cie-^''6D{xr-^6D{yr-^ < PD{t,x,y) < C2e-^i%(x)"-i5i5(y)-\ 

where — Ai < is the largest eigenvalue of the -generator of X. 



By integrating the above two-sided heat kernel estimates in Theorem 17.21 with respect to t, 
one can easily obtain the following sharp two-sided estimate of the Green function G£)(x,y) = 
Jo° Poit, x,y)dt of a censored stable process in a bounded C^'^ open set D. 
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Corollary 7.3 ([.16^ Corollary 1.2]) Suppose that n > I, a £ (1,2) and D is a bounded C^'^ 
open set in R". Then on D x D, we have 

a-1 

I J, ( 1 A ''oy-^>''oyy> ' 



Sharp two-sided estimates on the Green function are very important in understanding deep 
potential theoretic properties of Markov processes. When D is a bounded C^'^ connected open sets 
in R"" and n > 2, estimates in Corohary 17.31 had been obtained in [13j . 

Our approach in [16] is adapted from that of [15]. In [15j . the following domain monotonicity 
for the killed symmetric stable processes is used in a crucial way. Let Z he a symmetric a-stable 
process and be the subprocess of Z killed upon leaving an open set D. If U is an open subset 
of D, then is a subprocess of Z^ killed upon leaving U. However censored stable-like processes 
do not have this kind of domain monotonicity. So there are new challenges to overcome when 
studying heat kernel estimates for censored stable processes. A quantitative version of the intrinsic 
ultracontractivity, a crucial use of boundary Harnack inequality for censored stable process and the 
reflected stable process X all played an important role in our approach in [16] . 



8 Concluding Remarks 

In this paper, we surveyed some recent progress in the study of fine potential theoretic properties 
of various models of symmetric discontinuous Markov processes that the author is involved. To 
keep the exposition as transparent as possible, sometimes we did not state the results in its most 
general form. For example, results in Sections 2 and 3 hold for general d-sets F in R" and for 
F being a measure- metric space satisfying certain conditions, see [181 119): and the Dirichlet heat 
kernel estimates in Section 7 in fact holds also for censored stable-like processes, see [16]. Two-sided 
transition density function estimates for relativistic stable processes in C^'^ open sets have recently 
been established in [17]. The study of sharp two-sided heat kernel estimates for discontinuous 
Markov processes is in its early stage and is currently a very active research area. There are many 
questions waiting to be answered and several active studies are currently underway. For instance, 
it is natural to study the large time estimate for p{t, x, y) of the processes considered in Section 3 
for the case of 72 > 71 > 0. It is also nature and important to investigate the sharp two-sided heat 
kernel estimates for stable processes of mixed type in C^'^-open sets, and for Levy processes that is 
the independent sum of a Brownian motion and a symmetric stable process in C^'^-domains. Some 
promising progress has already been made in these studies. 



References 

[1] M. T. Barlow, R. F. Bass, Z.-Q. Chen and M. Kassmann, Non-local Dirichlet forms and 
symmetric jump processes. Trans. Amer. Math. Soc. 361 (2009), 1963-1999. 



23 



[2] M. T. Barlow, A. Grigor'yan and T. Kumagai, Heat kernel upper bounds for jump processes 
and the first exit time. J. Reine Angew. Math. 626 (2009), 135-157. 

[3] R. F. Bass, SDEs with jumps. Notes for Cornell Summer School 2007. 
|http:/ /www.math.uconn.edu/~bass/cornell.pdf| 

[4] R. F. Bass and D. A. Levin, Harnack inequalities for jump processes. Potential Anal. 17 (2002), 
375-388. 

[5] R. F. Bass and D. A. Levin, Transition probabilities for symmetric jump processes. Trans. 
Amer. Math. Soc. 354 (2002), 2933-2953. 

[6] I. Benjamini, Z.-Q. Chen and S. Rohde, Boundary trace of reflecting Brownian motions. Probab. 
Theory Relat. Fields, 129 (2004), 1-17. 

[7] J. Bertoin. Levy Processes. Cambridge University Press, 1996. 

[8] K. Bogdan, K. Burdzy and Z.-Q. Chen, Censored stable processes, Probab. Theory Related 
Fields, 127 (2003), 89-152. 

[9] R. M. Blumenthal and R. K. Getoor, Some theorems on stable processes. Trans. Amer. Math. 
Soc. 95 (1960), 263-273. 

[10] L. A. Caffarelli, S. Salsa and Luis Silvestre. Regularity estimates for the solution and the free 
boundary to the obstacle problem for the fractional Laplacian. Invent. Math. 171(1) (2008) 
425-461. 

[11] E. A. Carlen and S. Kusuoka and D. W. Stroock. Upper bounds for symmetric Markov 
transition functions. Ann. Inst. Heri. Poincare- Probab. Statist., 23 (1987), 245-287. 

[12] Z.-Q. Chen, Multidimensional symmetric stable processes. Korean J. Comput. Appl. Math. 6 
(1999), 227-266. 

[13] Z.-Q. Chen and P. Kim, Green function estimate for censored stable processes, Probab. Theory 
Related Fields, 124 (2002), 595-610. 

[14] Z.-Q. Chen, P. Kim and T. Kumagai, Weighted Poincare inequality and heat kernel estimates 
for finite range jump processes. Math. Annalen 342 (2008), 833-883. 

[15] Z.-Q. Chen, P. Kim and R. Song, Heat kernel estimates for Dirichlet fractional Laplacian. To 
appear in J. European Math. Soc. 

[16] Z.-Q. Chen, P. Kim and R. Song, Two-sided Heat Kernel Estimates for Censored Stable-like 
Processes. To appear in Probab. Theory Related Fields. 

[17] Z.-Q. Chen, P. Kim and R. Song, Relativistic Stable Process and its Transition Density Esti- 
mates in Unbounded Open Sets. Preprint, 2009. 

[18] Z.-Q. Chen and T. Kumagai, Heat kernel estimates for stable-like processes on d-sets. Stochastic 
Process. Appl, 108 (2003), 27-62. 

[19] Z.-Q. Chen and T. Kumagai, Heat kernel estimates for jump processes of mixed types on 
metric measure spaces. Probab. Theory Related Fields 140 (2008), 277-317. 

[20] Z.-Q. Chen and T. Kumagai, A priori Holder estimate, parabolic Harnack inequality and heat 
kernel estimates for diffusions with jumps. To appear in Revista Mathemdtica Iberoamericana. 

[21] Z.-Q. Chen and R. Song, Estimates on Green functions and Poisson kernels of symmetric 
stable processes. Math. Ann., 312 (1998), 465-601. 



24 



[22] Z.-Q. Chen and R. Song, Drift transforms and Green function estimates for discontinuous 
processes. J. Punct. Anal., 201 (2003), 262-281. 

[23] Z.-Q. Chen and J. Tokle, Sharp transition density estimates for symmetric stable processes in 
unbounded open sets. In preparation. 

[24] I. W. Herbst and D. Sloan, Perturbation of translation invariant positive preserving semigroups 
on L2(R^). Trans. Amer. Math. Soc, 236 (1978), 325-360. 

[25] S. R. Hurst, E. Platen and S. T. Rachev, Option pricing for a logstable asset price model. 
Math. Comput. Modelling 29 (1999), 105-119. 

[26] N. Ikcda, M. Nagasawa and S. Watanabc, A construction of Markov process by piecing out. 
Proc. of Japan Academy 42 (1966), 370-375. 

[27] A. Janicki and A. Weron, Simulation and Chaotic Behavior of a-Stable Processes. Dekker, 
1994. 

[28] J. Klafter, M. F. Shlesinger and G. Zumofen, Beyond Brownian motion. Physics Today, 49 
(1996), 33-39. 

[29] V. Kolokoltsov, Symmetric stable laws and stable-like jump-diffusions. Proc. London Math. 
Soc. 80 (2000), 725-768. 

[30] T. Kulczycki, Properties of Green function of symmetric stable processes, Probab. Math. Stat. 
17(1997), 381-406. 

[31] A. Matacz, Financial modeling and option theory with the truncated Levy process. Int. J. 
Theor. Appl. Finance 3(1) (2000), 143-160. 

[32] P.-A. Meyer, Renaissance, rccoUcments, melanges, ralentissement de processus de Markov. 
Ann. Inst. Fourier, 25 (1975), 464-497. 

[33] G. Samorodnitsky and M. S. Taqqu, Stable Non-Gaussian Random Processes. Chapman &i 
Hall, New York-London, 1994. 

[34] L. Silvestre, Holder estimates for solutions of integro differential equations like the fractional 
Laplace. Indiana Univ. Math. J. 55 (2006), 1155-1174. 

[35] R. Song and Z. Vondracek, Parabolic Harnack inequality for the mixture of Brownian motion 
and stable process. Tohoku Math. J., 59 (2007), 1-19. 

[36] D. W. Stroock, Diffusion semigroup corresponding to uniformly elliptic divergence form oper- 
ator. Lect. Notes Math. 1321, 316-347, Springer- Verlag, 1988. 

[37] Q. S. Zhang, The boundary behavior of heat kernels of Dirichlet Laplacians, J. Differential 
Equations, 182 (2002), 416-430. 



Zhen-Qing Chen 

Department of Mathematics, University of Washington, Seattle, WA 98195, USA and 
Department of Mathematics, Beijing Institute of Technology, Beijing 100081, P. R. China 
E-mail: zchen@math.washington.edu 



25 



